Testing the interaction between dark energy and dark matter via latest observations 
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Cosmological analysis based on currently available observations are unable to rule out a sizeable 
coupling between dark energy and dark matter. However, the signature of the coupling is not easy 
to grasp, since the coupling is degenerate with other cosmological parameters, such as the dark 
energy equation of state and the dark matter abundance. We discuss the possible ways to break 
such degeneracy. Based on the perturbation formalism, we carry out the global fitting by using 
latest observational data and get a tight constraint on the interaction between dark sectors. We 
find that the appropriate interaction can alleviate the coincidence problem. 
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PACS numbers: 98.80.Cq 



I. INTRODUCTION 



There has been convincing evidence indicating that 
our universe is composed of nearly 25% cold dark matter 
(DM) plus a small fraction of baryonic matter and around 
70% dark energy (DE) ^] . One leading candidate of such 
a DE is the cosmological constant, representing a vacuum 
energy density with constant equation of state w = —1. 
However, it is difficult to understand such a cosmological 
constant in terms of fundamental physics. Its observed 
value is far below that estimated in quantum field theory, 
what is referred to as the cosmological constant problem. 
Moreover, using the cosmological constant to explain the 
DE, there is no natural understanding why the constant 
vacuum energy and matter energy densities are precisely 
of the same order today. This is the so-called coincidence 
problem. 

Considering that DE and DM are dominant sources of 
the content of the universe, it is natural, in the frame- 
work of field theory, to consider the inevitable interaction 
between them [2|. An appropriate interaction between 
DE and DM can provide a mechanism to alleviate the 
coincidence problem Q-Q- A non- minimal coupling in 
dark sectors can affect significantly the expansion history 
of the universe and the density perturbation evolution, 
changing the growth history of cosmological structures. 
The possibility that DE and DM interact with each other 
has been widely discussed recently A number of 

studies have been devoted to analyze the constraints on 
the dark sectors mutual interaction from the probes of 
the cosmic expansion history by using the WMAP, SNIa, 
BAO and SDSS data etc [i3-[23|. Interestingly it was 
disclosed that the late ISW effect has the unique ability 
to provide insight into the coupling between dark sectors 

Furthermore, complementary probes of the coupling 
within the dark sectors have been carried out in the study 
of the growth of cosmic structure [l^l-fl^]- It was found 
that a non-zero interaction between dark sectors leaves 
a clear change in the growth index [2^ [2^. In addi- 



tion, it was suggested that the dynamical equilibrium 
of collapsed structures such as clusters would acquire a 
modification due to the coupling between DE and DM 
(26l [27j . Comparing the naive virial masses of a large 
sample of clusters with their masses estimated by X-ray 
and by weak lensing data, a small positive coupling has 
been tightly constrained [27], which agrees with the re- 
sults given in [19] from CMB. The small positive cou- 
pling indicates that there is energy transfer from DE to 
DM, which can help to alleviate the coincidence problem 

Mm- 

Both DE and DM are currently only detected via their 
gravitational effects and any change in the DE density is 
conventionally attributed to its equation of state w. This 
status leads to an inevitable degeneracy while extracting 
the signature of the interaction between dark sectors and 
other cosmological parameters. In this work, we will first 
discuss the degeneracy between the DE and DM coupling 
and the equation of state (EoS) of DE in the background 
dynamics. Furthermore, we will extend our discussion to 
the perturbed spacetime by considering the perturbation 
evolution of DE and DM. We review the formalism of the 
perturbation theory when there is an interaction between 
dark sectors. Based upon this formalism we explore the 
possibility of breaking the degeneracy between the cou- 
pling and other cosmological parameters, such as the EoS 
parameter w of DE as well as the DM abundance. This 
can help us extract a tighter constraint on the interaction 
between dark sectors from observations. 



II. THE BACK GROUND DYNAMICS 

In the spatially flat Friedmann-Robertson- 
Walker(FRW) background, if there is an interaction 
between DE and DM, neither of them can evolve 
independently. The (non) conservation equations are 
described by 



p'^ -I- 3^(1 -f w)pd 



aQa 



(1) 
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where the subscript "c" represents DM and "d" stands 
for DE. Q\ is the term leading to energy transfer. Con- 
sidering that there is only energy transfer between DE 
and DM, we have Qc — —Qd = Q- The sign of Q de- 
termines the direction of energy transfer. For positive 
Q, the energy flows from DE to DM. For negative Q, 
the energy flow is reversed. Since we know neither the 
physics of DM nor that of DE at the present moment, 
we cannot write out the precise form of the interaction 
between them from first principles (see, however, (H). 
One has to specify the interaction either from the outset 
4 1 , or determine it from phenomenological requirements 

aliil. 

For the sake of generality, we consider the phenomeno- 
logical description of the interaction between DE and DM 
as a linear combination of energy densities of dark sec- 
tors Q — 3H{^ipm + £,2Pd) 0, Ell- In studying the cur- 
vature perturbation it has been made clear that when 
the interaction is proportional to the energy density of 
DE {Q = 3H(^2Pd), we get a stable curvature perturba- 
tion; however, when the interaction is proportional to 
the DM density {Q = iH^ipm) or total dark sectors 
{Q = 3H^{pm-\-pd))^ the curvature perturbation can only 
be stable when the constant DE EoS satisfies i/; < — 1 Q- 
For the case of a time-dependent DE EoS, the stability 
of curvature perturbations was discussed in [loj . 

The presence of the coupling also changes the DM and 
DE redshift dependence acting as an extra contribution 
to their effective equation of state. Indeed, the effective 
background equations of state for the two fluids are 



Wc,eff 



?^d,cff 



aQc 
3'Hpc 
aQd 



SHpd 



(2) 



Choosing different forms for the interaction, the effective 
background EoS can be read from table I. For conve- 
nience, we label our models with Roman numbers. 

We define r = Pc/pd as the ratio of the energy densities 
of DM and DE. In order to solve the coincidence problem 
we require the ratio of r to be a constant in the expansion 
history of our universe. This leads to the condition r' = 

Pd ^ ''pd ^ which yields a quadratic equation. 



(3) 



When the coupling is proportional to the dark energy 
density (Models I, II) — 0, Eq. ^ has only one root, 
— — ' which will appear in the future. In 21] it 
was found that when the coupling is proportional to the 
energy density of DE, ^2 > — 2K;r2c leads to a negative 
energy density of cold DM(pc < 0) in the past, which 
is unphysical. We thus require ^2 < —2wflc- In such a 
case, the effective DE EoS reads Wd,cff = w + £,2, which 
shows that the coupling and the DE EoS are entangled 
and it is impossible to distinguish the coupling from DE 
EoS by investigating the background evolution. 



When the coupling is proportional to the dark matter 
density or to the total dark sector (Models III, IV), there 
are two roots of the quadratic equation Eq. 



1 



1 



Ki)i = -^{w + a + 6) + 2 v/(a + 6 + - 466 



Hi)2 = -^{w + 6) - 2 v/(a + 6 + wr - 466 

The first root happens in the past and the second one 
happens in the future. For the interaction proportional 
to the dark matter density (Model III) (6 = 0), when 
6 < —w/A, the roots are real and can be simplified. We 
have, 



(r6)i = -(i« + 6) >o 
Hi)2 - 



(4) 



The effective DE EoS can behave differently in the past 
and in the future. In the early time of the universe, the 
effective DE EoS can be of the form 



Wd,cS 



<1 



(5) 



We learn that the coupling is entangled with the DE EoS. 
When the universe evolves to the present time or into the 
future, the effective DE EoS reads 



w 



(6) 



which does not depend on the coupling and thus the de- 
generacy between DE EoS and coupling no longer exists 
in the background dynamics. 

When the interaction is proportional to the energy den- 
sity of total dark sectors (Model IV), (6 =6), Eq- ® 
has two real roots when ^ < —w/A, 



H)2 



1 



1 



-iw + 20 



(7) 



according to which the effective DE EoS turns out to 
be Wd,cfi ~ in the early time of the universe and 
Wd,cS ~ w + £,'^ + at late times. There is a degener- 
acy between the coupling between dark sectors and the 
DE EoS in the background dynamics. From the back- 
ground dynamics we see that when we introduce the in- 
teraction between DE and DM, it is possible to have the 
scaling solution of the ratio between DM and DE, which 
can help to alleviate the coincidence problem. However, 
in the background dynamics there appears an inevitable 
degeneracy between the coupling between dark sectors 
and the DE EoS. In general this degeneracy cannot be 
broken by just investigating the dynamics of the back- 
ground spacetime, except in the case when the coupling 
is proportional to the dark matter density (Model III). 
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TABLE I: Coupling models 



Model 


Q 


DE EoS 


U)c,eff 


lt)d,cft 


Constrains 


I 


3C2i?Pd(Ci = 0) 


-1 < -u; < 




W+^2 


^2 < -2wQ.c 


II 


3C2i?Pd(Ci = 0) 


Ul < — 1 




W+C2 


^2 < -2wQ.c 


III 


HiHp,{i2 = 0) 


w < —I 




w + 5ir 


< ^1 < -w/A 


IV 




w < —I 






< C < -1^/4 



In the following we are going to explore the possibility of 
breaking the degeneracy between the coupling and other 
cosmological parameters in the perturbed spacetime by 
considering the perturbation evolution of DE and DM. 

III. THE PERTURBATION FORMALISM 

In this section, we will go over the first order met- 
ric perturbation theory in the presence of coupling be- 
tween DE and DM. The perturbed space-time at first 
order reads 

+ ( f + 20) (5,y- dx' dx^ + D,j Edx' dx'], (8) 

where ip^B^cji^E represent the scalar metric perturba- 
tions, a is the cosmic scale factor and 

A, = {d,d, - i%V2). (9) 



We work with the energy-momentum tensor 

T^"" = pV^U" + pig'"' + WU"), (10) 

for a two-component system consisting of DE and 
DM. The covariant description of the energy-momentum 
transfer between DE and DM is given by [28| 

V^T(,)'^- = Q(,)'' (11) 

where Q(a)'' is a four vector governing the energy- 
momentum transfer between different components. The 
subindex A refers to DM and DE respectively. For the 
whole system, the energy and momentum are conserved, 
and the transfer vector satisfies 

A 

In the Fourier space the covariant form of perturbed 
Eq. (HB) reads [1 ^ 



opx 



(vx + BY + nil-3wx){vx + B) - 



k 5px 
1 + WA Spx 



-3{l + wx)ck' + {2i,-Sx)—^ + 



Sx = -- 



1 + WA 



(wA + B) + k^- 



a 



2^)0 
A 

PA 



PA 



wxa'^Ql 



wx)px 



■B + 



(12) 

a'^SQpx 
(1 -I- wx)px 
(13) 



By constructing gauge invariant quantities [Ji 



SQix^SQpx-Ql^; 



2fc' 



Ql 



1 



6 ' 



Dgx = Sx- 



Pa 
pxH 



E' 

^^="^"2fc' 



(14) 



we obtain the general gauge-invariant perturbation equa- 
tions for DM and DE respectively. 



u' + nUc 



Pc 



(15) 




- iw' + iicl - w)^ + -^7^ I $ + \m{cl -w) + 

Pd PdH pd j L 



Pd 



= -kUd + m{ct-ci 



2^Pd 



Pd (1 + w)k 



2* 



01 



a'Q'd^ -'Q'df^ 
PdT-L Pd Kt-L 



Pd 



D 



gd 



PdT-L 



U'd + n{l - 3w)Ud - kC',D,d + kC!^<i> ~ {Cl - Cl) + (1 + w)k^ - '^Ud + {l + w) 

Pdti ' \ + w Pd Pd 



pd 



Pd 



(16) 



where Dgx is the energy density contrast in spatial flat obtain the perturbation form 



gauge[28], U\ = {l + w\)V\ and V\ is the gauge invariant 
peculiar velocity. We have employed 



Spd 

Pd 



Cl5d-{Cl-Cl)'- 



■ Vd + B 



Pd 



(17) 



where is the effective sound speed of DE at its rest 
frame which is gauge invariant under gauge transforma- 
tion and Cl is the adiabatic sound speed. 

Eqs. ([T5)) and p6|) are the most generic form regardless 
of the detailed description of the interaction Q(a)''. 

The four vector Q(a)'^ can be phenomenologically de- 
composed into two parts with respect to a given observer 
A' with four velocity C^(a')^- 



>,(\) 



}(\\')U(yf 



(18) 



where (3(aa') = ^U(\i)^Q(\)'^ is the energy transfer rate 
of A component observed by A' observer. i^(AA')^ — 
h[\i)^ ^Q(\)^ is the corresponding momentum transfer. 
In appendix IVI Bl we show that such decomposition of 
Q{\)' and its perturbed form are identities regardless of 
the observer. The four vector Q(\)^ must be specified di- 
rectly, according to the physical meaning. As discussed 
in IVI A| in cosmology, we need to specify the coupling 
vector Q(\)^ in the co-moving frame as 



<3(A)' 



Q{\) 



,0,0,0 



(19) 



Q{X) is the module of four vector Q[x)^ . The perturbed 
form <5Q?x-, can be uniquely determined from the back- 



'(A) 

ground energy- momentum transfer Q^x) • From 

Q(\fQ(X)ti = 50o('9(A)")^ = ^Q(A)^ 



(20) 



where Qja) = aQ(A)'^ is a scalar in the FRW space and 
the minus sign indicates that Q{\)^ is time-like, we can 



(21) 



It can be shown that the zero component of the per- 
turbed energy-momentum transfer is covariant. The 
spatial component of the perturbed energy-momentum 
transfer SQ^^^^^ is independent of the zeroth component. 
It refers to the non-gravitational force and is composed 
of two parts, 

SQpx - SQ'px It +Q?A)"* (22) 

where SQp\ is the potential of the perturbed energy- 
momentum transfer SQ^^^y SQpx It external non- 
gravitational force density and Vt is the average velocity 
of the energy transfer, vt is a free quantity which needs to 
be specified according to physics. In Q, was allowed to 
follow the peculiar velocity of DM or DE respectively. If 
we allow Vt = Vc or vt — Vd, we can reproduce the result 
in 8]. In our analysis, we consider the fact that there 
is no non-gravitational interaction in the DE and DM 
coupled system, only the inertial drag effect appears in 
the system due to the stationary energy transfer between 
DE and DM as discussed in 13). Thus we set Vt — and 
SQpx |t= Oi which leads to the vanish of JQ'- 

In constructing the four vector Eq. (jl9[) . the mod- 
ule Qx can be chosen as any combinations of scalar 
in the FRW space, such as the energy density p(A) — 
T(\fU(^x)iJJ(\)V, expansion i7(A) = \'^t,U(xf, or any 
scalar function thereof. Considering that Qx is observer 
independent and so does the energy density and its per- 
turbed form, we require to be a global quantity to 
avoid the ambiguity of the observer dependence. In gen- 
eral phenomenological description, we can assume the 
coupling to be given by 

Qc = Q = 3H{^ip,+^2Pd) 

Qd = -Q = -3i7(aPc + 6Pd) • 
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The perturbed forms read 

a a 

a a 
I 

Pc 

Pd 



where r — pd Pd- Inserting the above expressions into 
Eqs. (IT5|) . ([TE)) and neglecting the spatial perturbations 

I 



The general gauge invariant formalism fully removes the 
ambiguity of gauge choice. However, it can be solved 
through gauge-dependent methods by picking a peculiar 
gauge, without loosing generality (see.[23| chapter III for 
details). The results will be the same for different gauges 
if the gauge is fully fixed (see appendix IVI CI) . In the 
following discussion we will choose the Conformal Newto- 
nian gauge together with the adiabatic initial conditions 
as we used in [9]. 

Now we are in a position to use the perturbation for- 
malism to study the influence of the interaction between 
dark sectors and other cosmological parameters on the 
CMB power spectrum. In Figs[T]and[2]we illustrate the 
theoretical computation results of the CMB power spec- 
trum for different cosmological parameters. 

Fixing the DM abundance, we see in the CMB TT 
angular power spectrum (FigslT]) that the change of the 
constant DE EoS only modifies the low-1 part of the spec- 



The first terms in the last two equations were omitted in 
our previous work fl9l |. however, only slightly modifica- 
tion will be brought by adding these terms. 

The gauge invariant quantities 5Q°J and 8Q^J ap- 
peared in Eqs. P3|) . as defined in [l^ can be cal- 
culated as. 



^Qp\ = 0' obtain the general gauge-invariant pertur- 
bation equations for DM and DE respectively. 



(23) 



(24) 



trum while leaves the acoustic peaks almost unchanged. 
When the constant DE EoS w > —1, the low-1 spec- 
trum gets enhanced with the increase of the value of w. 
Such a property keeps valid when the DE EoS parame- 
ter is a constant smaller than -1, namely the phantom 
case w < —I. However when the DE is in the phan- 
tom region, the enhancement of the low-1 spectrum due 
to the increase of the w is less sensitive than that of the 
quintessence DE. 

In the low-1 CMB spectrum, we see from FigsH] that 
the coupling between dark sectors can also change the 
spectrum. As the coupling becomes more positive, the 
low-1 spectrum is further suppressed. When the inter- 
action between dark sectors is proportional to the dark 
matter or total dark sector energy density, the low-1 spec- 
trum is more sensitive to the change of the coupling than 
the DE EoS. 



~3H(a + 6/0* + iUiCiDg, + 6^gd/r}$ + 3(6- 



6 Pd^ 

r Pd' 



Pc H Pc 



n 



= 3H(6'^ + 6)«'-3H{6^, 



gel 



6i^s4 - 3(6^-^ + -6)* - 

Pc Pd 



P'c 1 P'dc^\^ ^ I " 



Pd 



Pd 



n 



D' 



gc 
U' 



D' 



gd 



= -kUc + 3H*(6 + 6/0 - 3(6 + 6/0*' + SnUDgd - Dgc)/r 

- -nUc + k^-3n{^i+^2/r)Uc , 

- -m{C^ - w)Dgd + {3w' - 9H{w - Cl) (6r + 6 + 1 + ^f)} $ 

- m^{Cl - Cl)^ + 3(6r + 6)*' - 3*H(6r + 6) + 3n^ir{Dgd - D 



U'd 



gc) 



9'H\Cl 



Ud 



(1 + w)k 



-nil - 3w)Ud - 3kCl (6r + 6 + 1 + w) * + 3H(Cf - Cl){^ir + 6) 

+ i{cl - cDnUd + kclDgd + (1 + w)fc* + 3H(6r + 6)c/d- 

I 



Ud 
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FIG. 1: The CMB TT power spectrum for different interaction models between DE and DM. 
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FIG. 2: The dependence of CMB angular power spectrum on physical cosmological parameter = ficft^ 



Beside the low-1 CMB spectrum, the interaction be- 
tween DE and DM can also influence the acoustic peaks. 
This feature is interesting, since this property differs from 
that of the DE EoS and can help to break the degeneracy 
between the interaction between dark sectors and the DE 
EoS. 

The above discussion is valid for fixed DM abun- 
dance. Now we investigate the dependence of CMB an- 
gular power spectrum on the abundance of cold DM, 
Wc = ilch^. Although the abundance of the DM does 
not affect much on the low-1 CMB power spectrum, it 
quite influences the amplitude of the first and second 
acoustic peaks in CMB TT angular power spectrum (see 



Fig. [2]). Decreasing the cold DM abundance lUc will en- 
hance the acoustic peaks. This effect is degenerated with 
the influence given by the dark sectors' interaction as we 
observed in FiglTJ A possible way to break this degener- 
acy is to consider the influence of the interaction on the 
low-1 CMB spectrum. Moreover, we can include further 
observations to get a complementary constraint on the 
DM abundance and this in turn can help to constrain 
the coupling between dark sectors. 

In order to extract the signature of the interaction and 
constraints on other cosmological parameters, we need 
to use the latest CMB data together with other observa- 
tional data. We perform such a task in the next section. 
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IV. GLOBAL FITTING AND COSMOLOGICAL 
COINCIDENCE PROBLEM 

In this section we confront our models with observa- 
tional data by implementing joint likelihood analysis. We 
take the parameter space as 

P = (/i,tJ6,iUcdm,T,ln[10^°A5],ns, 6,6,1") 

where h is the bubble constant, = il^h'^ ,ujcdm — 
^cdmh^ , As is the amplitude of the primordial curvature 
perturbation, Us is the scalar spectral index, 6 and 6 are 
coupling constants proportional to the energy density of 
DM and DE respectively, w is the DE EoS. We choose the 
flat universe with fli^ — and our work is based on CM- 
BEASY codeflil. We use the CMB anisotropy data from 
the seven-year Wilkinson Microwave Anisotropy Probe 
(WMAP). The Priors for Baycsian Analyses are pre- 
sented in table HIl The fitting results in the la range 
are listed in table IIIII We plot the likelihood for ilch^, 
the coupling between dark sectors and the DE EoS in Fig 
|3l It is clear that, when the interaction is proportional 
to the energy density of DE, CMB data alone can not 
impose good constraints on ilch^, the coupling and the 
DE EoS all together. This can be explained by our theo- 
retical analysis shown in FigslT]and[5] At the low-1 CMB 
spectrum, it is impossible to distinguish the DE EoS and 
the coupling influences; in acoustic peaks it is hard to 
break the degeneracy between the coupling and the DM 
abundance. 

When the interaction between dark sectors is propor- 
tional to the energy density of DM or total dark sectors, 
CMB data alone can impose tight constraints on cou- 
plings and flch^, but it can not impose good constraint 
on the DE EoS w. This result can also be understood 
from our analysis in FigsH] and [D The degeneracy be- 
tween the coupling and the DM abundance can be broken 
by looking at the low-1 CMB spectrum, since the CMB 
spectrum is not sensitive to the change of the DE EoS 
in the low-1 spectrum. In order to get tighter constraint 
on D.chF' , we use the BAO distance measurements [sO] 
which are obtained from analyzing clusters of galaxies 
and tests a different region in the sky as compared to 
CMB. BAO measurements provide a robust constraint 
on the distance ratio 



rs{zd)/Dy{z) 



(25) 



where D^{z) = [{1 + z)^D\z/H{z)]'^/Hs the effective dis- 
tance Da is the angular diameter distance, and H{z) 
is the Hubble parameter. rs{zd) is the comoving sound 
horizon at the baryon drag epoch where the baryons de- 
coupled from photons. We numerically find Zd using the 



condition t/R = 1,R 
Xbao calculated as 30], 



3 Pb. 
4 



as defined in 



The 



o - -^obs ^ rr' 1 , ->obs , 

x|^o = (d-d fC-\d-d ) (26) 
where d = (^=0.2,^^=0.35)^, d"*" = (0.1905,0.1097)^ 



and the inverse of covariance matrix (3C 



30124 -17227 
-17227 86977 



Furthermore, we add the BAO A parameter [3E 



(27) 



i;(0.35)i/3 



0.35 



0.35 



dz 

m 



2/3 



0.469(ns/0.98)-''-^^± 0.017 



(28) 



where E{z) = -^^^ and are the scalar spectral index. 
In order to improve the constraints on the DE EoS w, 
we use the compilation of 397 Constitution samples from 
supernovae survey [s^l- We compute 



E 



(29) 



and marginalize the nuisance parameter. In addition to 
the above mentioned data sets, we also add the latest 
constraint on the present-day Hubble constant [35] , 



Ho = 74.2 ± 3.6kms"^Mpc"\ 



We implement the joint likelihood analysis that is. 



(30) 



= XwMAP + XsN + Xbao + Xh,, ■ (31) 



The fitting results are shown in Table HVl The cosmolog- 
ical parameters are well constrained. When the coupling 
between dark sectors is proportional to the energy den- 
sity of DE, its value is constrained up to a few percent. 
When the coupling is proportional to the energy density 
of DM or total dark sectors, its constraint is pretty good 
and reads 6 = 0.0006 



+0.0006 J c 
-0.0005 ? 



0.0006 



-1-0.0005 
-0.0006 



re- 



spectively. In 1 CT range the couplings are positive. The 
likelihoods of the fitting results for the DM abundance, 
DE EoS and the coupling between dark sectors are shown 
in Fig m Compared with the WMAP data alone, we see 
that the joint analysis by including other observational 
data provides tighter constraints on the cosmological pa- 
rameters. 

The positive coupling can help us to alleviate the cos- 
mological coincidence problem (Tol . [2l| . As shown in Fig 21 
the energy density of DE and DM in standard ACDM 
model are only comparable at present moment. The thick 
black line representing the quantity logiQ^pc/ Pd) is lin- 
early proportional to log^Q a and almost precisely crosses 
the origin in the whole expansion history of the universe. 
However, this is hard to be convincing and achieving. If 
we want the energy density of DM to be comparable to 
that of the DE at the present moment, such an origin 
crossing can only be realized by tuning the initial con- 
ditions at the early time of the universe over 30 orders 
in energy density contrast r — pd Pd- If there is little 
change in the initial condition, the log]^g(pc/pd) cannot 
cross the origin and at the present the energy densities 
of DM and DE cannot be comparable. This problem can 
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be overcome by introducing the interaction between DE 
and DM. As an example we show the model when the 
interaction between DE and DM is proportional to the 
energy density of the total dark sectors (Model IV). In 
this case there are two attractor solutions of the ratio r 
during the expansion history of the universe 



1 



(32) 



by considering ^ is a small value from the fitting results, 
ri happened in the past and r2 will occur in the future. 
The behavior of the attractor solutions of the ratio r only 
depends on the coupling constant ^ and does not depend 
on the initial conditions at the early time of the universe. 
To see this point more clearly, we show that in Fig [SI the 
purple lines represent the density evolution of cosmolog- 
ical model with different initial conditions. The density 
contrast r at present is different for different initial con- 
ditions. However, they are all bounded in two attractor 
solutions ri ~ ^, r2 ~ 1/^ in the plane. Adopting the 
coupling constant value from the fitting, 1/^ ~ 10"* and 
^ ^ 10"'', we have the ratio r in the range 10"'' < r < 10"* 
during the universe history. Thus the change of the ratio 
r is much smaller than that of the ACDM model so that 
the period when the DE and DM are comparable is much 
longer than that of the ACDM model. The cosmological 
coincidence problem can thus be greatly alleviated. 
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VI. APPENDIX: 

A. covariant equation of motion for interacting 
system 

The basic dynamics of interacting systems in classical 
mechanics is described by Meshchersky's equation. 



dm , 



dt 



dt 



pa 



(33) 



where ^ is the rest mass (energy) transfer rate by the 
moving system, F"^ is the external force, is the 

momentum transfer and ^ttw" is the inertial force. 

at 

The classical Meshchersky's equation can be extended 
to special relativity and general relativity. In the frame- 
work of special relativity, it was first derived by Ackert 
[sgI ] and then summarized by Seifert in [sl] . Here we first 
extend their results on the equations of motion to a co- 
variant form, then generalize them to curved spacetime. 
We consider a moving system with rest energy density p 
and the energy-momentum tensor. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have reviewed the formalism of the 
perturbation theory when there is an interaction between 
DE and DM. We have proposed a way to construct the 
coupling vector in a self consistent manner both in the 
perturbed form and in the background. Based upon the 
perturbation formalism we have studied the signature of 
the interaction between dark sectors from CMB angular 
power spectrum. Theoretically we found that there are 
possible ways to break the degeneracy between the inter- 
action, DE EoS and DM abundance. This can help to 
get tight constraint on the interaction between DE and 
DM. 

We have performed the global fitting by using the CMB 
power spectrum data from WMAP7Y results together 
with latest SNIa, BAO and Hq data to constrain the in- 
teraction between DE and DM. When the interaction 
between DE and DM takes the form proportional to the 
energy density of DM and the total dark sectors, in la 
range the coupling is found to be positive. The tight pos- 
itive coupling indicates that there is energy flow from DE 
to DM, which can help us to alleviate the cosmological 
coincidence problem. 

The question of how to improve the model is now much 
related to find a field theory based model for the interac- 
tion and how to relate the model to the standard model 
of particle interactions. This is currently under study. 



T 



ah 



arjb 



(34) 



where C/° = (^)'' is the four velocity. The covariant 
form of the equation of motion can be given by. 



(35) 



where is the energy transfer four velocity and Q'^ is 
the external four force density acting on the system. For 
a given inertial observer Z° = (^)° in Minkowski space- 
time, the ordinary derivative operator db vanishes on 
dbZ'^ = 0. The time- like part of the above equation reads, 

ZadbT"'' = pZaU^+pZaA^ + ZaU^pO 



= -^(7P)-7P0 



(36) 



where A" — U^dbV^ is the four acceleration, dot denotes 
C/'Sfc = |:,e = aaf/^ 7 = -^a^/^7t = -ZaUf are 
Lorcntz-boost factors and ZaQp represents the energy 
transfer density observed by Z". 
The space-like part reads. 



= -^hpv") + p^v^Q 
dr 



(37) 
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Model I 








-ii^' 






1^ 

< < < < ^< , 


1 " 

<lV. < < < I < , < 




0.05 0.1 


0.15 -1 


-0.75-0.5-0.25 


-0.2 ( 


] 0.2 


Model II 


1 j/< 






1 1 




M 




0.05 0.1 


0.15 -2.5 -2 -1.5 -1 -0.2 0.2 


Model III 


i 




1 1 


\ 


_'' 1 _ 




0.05 0.1 


0.15-2.5 -2 -1.5 -1 


0.005 0.01 


Model IV 


J 


1 




1 

-J 






' 1 _ 





0.05 p.1 0.15 -2.5 -2 -1.5 -1 0.005 0.01 

n h- w \ 



FIG. 3; The likelihood of cold dark matter abundance flch^, dark energy EoS w and couplings ^ for different cosmological 
models. The black lines denote the results obtained from WMAP 7 year data set alone and the red lines denote the results 
obtained from the combination of WMAP + SN + BAG + Ho 



TABLE II: Priors for Bayesian Analyses 

< Qch'' < 0.5 
-l<w< -0.1(Model I) -2.5 < w < -l(Model II,III,IV) 
-0.4 < 5 < 0.4 (Model I,II) < ^ < 0.02(Model,III,IV) 



where ha'^ — Sa'^ + ZaZ'^ is the projection operator, v'^ = 
ha'^U"' /"f,v^ = haWl l^t are three velocities observed by 
and ha'^Qp is the three force density acting on the 
moving system. If there is no expansion = in the 
system, Egs. ([361) and dST]) go back to Eqs.(125) and 
(126) in 

With the help of the covariant form, Eq. ([55]) can be 
generalized to curved spacetime by the "minimal substi- 
tution" db Vfc, 



(38) 



The above equation is the generalized Meshchersky's 
equation which is the basic equation of motion for in- 
teracting systems in curved spacetime. 

In order to give a clear physical interpretation on this 
equation, we study the dynamics in terms of the distri- 
bution function. The energy-momentum tensor can be 



written as 37 1 



T 



'Cib 



(39) 



where / — f{x°',p°') is a distribution function, m is the 
rest mass for moving particles, — mW^ is the four 
momentum and — - is the Lorentz-invariant volume ele- 
ment on the positive-energy mass shell. Eq. (j38p can be 
presented in the form 



ab 



C[f]p' 



df* ad^P 
m——p . 

as m 



(40) 



If contracted with the given four velocity, the LHS of the 
above equation gives rise to the monopole and dipole of 
Boltzmann equations. The coupling vector on the 
RHS consists of two terms but with distinct physical 
meanings. The first term is the four force density pro- 
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TABLE III: The fitting results using WMAP 7-year data alone 



Parameters 


Model I 


Model II 


Model III 


Model IV 


h 




i.uy_o 26 


o.soiU- 


0.83+^;f3 


0.1,11^ 


o.o224±J!;;;;- 




o.o2i9tn;!i!^ 


0219+"-'*""'* 

U.UZ13_Q 0006 




< 0.111(68%CL) 


< 0.151(68%CL) 


n 1 1 7+u.u()y 

^•^^ ' -0.007 


o.ii9iii:r 


T 


0.084-'-- 


O.OSSiiJ;-^ 


0.087ini? 


o.o85tn;^ 


Us 


966+" "^* 


0.957t-,il 


"•^^^-0.016 


Q43+»-Ui''' 

U.9y:0_o 018 


ln[10'"As] 


O.U/ l_o.036 


O.U(Z_o.o35 




■J-"' ' -0.035 


w 


< -0.694(68%CL) 


unconstrained 


unconstrained 


unconstrained 








0010+"-""^^ 

U.UUiU_o 0010 


0011+"-""^" 

U.UUii_o.o011 



TABLE IV: The fitting results using WMAP + SN + BAO + Ho 



Parameters 


Model I 


Model II 


Model III 


Model IV 


h 


o.699±;!:;i- 


o.709t;!:;i- 


o.7oo±;!;;i;^3 


o.699t;i:[!i^3 




0.0224i-S 


n |-,ooo+U.UUU5 
U.UZZZ_o 0005 


U.UZZZ_o.ooo6 


U.UZZZ_o 0006 




o.io7t;!:r 




r, 1 1 n + U.UU3 
U.110_Q 


o.ii4t;!;r3 


T 


0.086±ni|: 


0.083±l!;;!i^ 


0.087±ni5 


0.087±;!;i]i^ 


Us 


u.yo/ _o 013 


U.yDl_0.013 


0.956±ro;i 




ln[10'"As] 


3 070+"-"^'' 


O.UU3_o 033 


o 074+0.038 


O.U/^_0.034 


w 


< -0.938(68%CL) 


1 r,q + U.IJ3 

l.UO_o 04 




_1 Qn + U.U3 
^^•'^'^-0.05 




-0.003±;!;i]^I 


024+" "^'* 

U.UZ'i_o 027 


0.0006±i!;r;^ 


0006+"-"'"'''' 

U.UUUD_o.o006 




FIG. 4: Cosmological coincidence problem 
duced by collisions, which is a space-like vector. 



Qf = / C[f]p 

J III- 

where C[f] is the collision kernel. The second term is the 
energy momentum transfer density along the direction of 
the average four velocity — QU^ is a time-like 
vector. 

QUt = / m^p^^ (41) 
J as m 

^ refers to the change rate in distribution function due 



FIG. 5: pc is the critical energy density today. The attractor 
solutions of r does not depend on the initial conditions at 
the early time of the universe. The purple lines represent the 
density evolution of cosmological model with different initial 
conditions. Noted by points, the density contrast r today are 
different for different initial conditions but they are bounded 
in two attractor solutions ri ~ ^, r-2 ~ 1/f in pc — Pd plane. 



to the varying rest mass of particles or varying comov- 
ing particle number density in the system investigated. 
However, it needs to be specified according to physics. If 
decomposed relative to a given observer with four veloc- 
ity Z", p"- can be represented as 



P 



m-i{p){Z''+v%p)) 



(42) 
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where ^{p) — —paZ"'/ni. Therefore, the thne-hke part of 
Eq. (gl]) reads 



(43) 



where 7 — —U^Za and 7(3 represent the energy transfer 
rate observed by The spatial part reads 



where is the average energy transfer velocity, 
^ j rn^-f{p)v''{p)(fip 



(44) 



(45) 



In hydrodynamics, ^Qvf refers to the viscosity due to 
the momentum transfer in different components. 

The above equations are quite general and now we con- 
centrate our discussion on the DE and DM coupled sys- 
tem. The external non gravitational force density acting 
on the system vanishes = in the background due to 
the homogeneous universe. Furthermore, noting that the 
spacetime is isotropic, ^ only depends on time and the 
contribution of v"" in Eq. psj) is counteracted in oppo- 
site directions, vanishes in the background. Only the 
energy transfer Q can be observed in the background. 

In the perturbed universe, as we have neglected the 
anisotropic stress-tensor, we assume that the perturba- 
tion is still isotropic, from Eq. (l45t we find that vf 
vanishes and the coupling vector is independent of the 
bulk motion of the component in the universe. 



B. covariant coupling vector in perturbed 
spacetime 

The general coupling vector Q[\)^ is independent of 
the choice of observer. But in the literatures, it is usually 
decomposed in two parts with respect to a given observer 
A' with four velocity U(x')'^- 

= -{Uix')^Qix)nU^X'f + hx')^^Qixf (46) 

where Q(aa') = ^U(\')^Q(X)'^ is the energy transfer 
rate of A component observed by the A' observer and 
F{xx')^ = ^(A')'^,y'3(A)'' is the momentum transfer ob- 
served by the A' observer, correspondingly. 

We can show that such a decomposition cannot bring 
substantial physics because Eq. (|46)) is an identity. Fur- 
thermore, we can show that the perturbed forms are also 
identities. The perturbation of the zero-th component on 
the RHS of Eq. reads 

{SQixx'))Uix')° + {Q(xx'))3U^X')° + SF(xx')'' 
= -myMx)''U(y)' - {U(y)JQix)nU(y)' 

-{Uix')^Qwnsu^x')' + shx')\Qw' 

+V)VQ(a)'='5Q(a)" (47) 



where we have used f7o(A') = "^i ~ ^/^^ ^^a')*" 



^Vf/(A')°, ^^A')0 



A')0 



in the derivation. We find 



that the nonzero SF(^xx')^ plays an important role in get- 
ting the identity in zero-th component. 

Similarly, the perturbation of the spatial component in 
the RHS of Eq. reads 



(^Q(AA'))t^(A')' + Q(AA')^C^(A')' + ^^(AA')' 



X')^Q(X)lU^y) 



+SihiX')\)Q'(x) + h^) 
= -(C/(AOoQ(A)°)<5C/(AO^ 
+h(x')\SQ(x)'' = SQ(x)' 



- iU(x'),' 

+ (C^(A')0< 



A') 



}x')SU^y)' 



(48) 



where C/(a')' = 0, (3(a')'' = 0, (5t/(A')' = d^V(y)/a have 
been employed. The net effect is /i(a')\<5Q(a)^ = SQix)' ■ 
Thus the i-th component is also an identity. 

Since Q{X)^ does not depend on observer, we need to 
specify such a coupling vector directly as discussed in 
IVI Al Once it is specified in the background, the zero-th 
component of the perturbed form SQ'^f^^) ''^^ uniquely 
determined by the background Q^^)- ^'^^ ^^^^ purpose, 
we consider the module of Q^x) 

Q(A)^Q(A)M = .900(Q(A)°)' = -Qix)\ (49) 

where Q{x) — o,Q{x)^ is a scalar on FRW space and the 
minus sign comes here because Q{x)^ is time-like. By 
considering the perturbation of the above equation, we 
find 



1 



Q(x) + -5Qix) = 5Ql 



(A) 



(50) 



where -0 arises from Sg^o. The first term is from the 
perturbation Scq = —ili/a, and the second term comes 
from the perturbation of the module. 

Q(x) is a scalar and under gauge transformation 



5Q{x) 
Noting [H 

we find 



JQ(A) - L^oQ^x) = 5Q{X) - Q[x)f- (51) 



(52) 



(A) 



(A) 



Q(A) 



(A) 



■SQ{x) 



??A)e"' 



(A) 



0' 



(53) 



which is consistent with the gauge transformation of 
(5(5(A) required by a covariant vector Q(a)'^[22]. 

The spatial part (5Q(a) independent of the zero-th 
component. It refers to non-gravitational force. The co- 
variant perturbation of the potential 5Qpx of the spacial 
part SQ'j^^^ can be written as, 



SQpx = SQ'px Ia' +Q(a)^^(a') 



(54) 
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where SQ^^^ \y is the perturbation SQp\ observed in A' 
rest frame and '9(>,)f(A') represents doppler effect, SQp\ 
satisfies the gauge transformation 



SQpx - 6Qpx + Q"(x)l3' 



(55) 



In DE and DM coupled system, we assume that SQ^^ \\' 
vanishes in the background SQp^ \back= 0, where 
V{back) = 0, and 



}p\ — SQpX \back +Q"\)V(back) = 



(56) 



There are no non-gravitational force and doppler effect 
produced by energy transfer. 



C. gauge conditions 

1. Conforrnal Newtonian gauge 

The conformal newtonian gauge is a set of coordinates 
in which the perturbed line element satisfies 



B = E = E' = 



(57) 



The gauge is fully fixed and thus the Barddeen's potential 
[39| can be simply calculated as 



The gauge condition fixes the expressions for the gauge 
transformation 



B 

T 



E' 
2fc2 



E 
2k 



(58) 



In particular, when we take the gauge transformation in 
different Conformal Newtonian coordinates B ~ B ^ 
E = E = E' = E' = 0, = P = 0, all the perturbations 
will have the same value eg. S(^x'j = 5(a), so that the 
Conformal Newtonian gauge yields unambiguous results. 



2. Gauge mode in Synchronous gauge and weak equivalence 
principle 

Synchronous gauge is defined by iJj = B — 0. The 
gauge invariant Barddeen's potential [391 in Synchronous 
gauge can be calculated by. 



(59) 



where rjx, are Synchronous gauge parameters, 

E^ 



VT 



(60) 



In contrast to Conformal Newtonian gauge, the met- 
ric conditions do not fully specify the gauge and need 
to be supplemented by additional definitions. When tak- 
ing the gauge transformation in Synchronous coordinates 
ip — tp — B = B = 0, it defines the gauge transformation 
up to two arbitrary constants Ci, C2. These constants 
manifest themselves in time and spatial coordinate trans- 
formation 



Cia 



f5 ^ ^kC, I - 
a 



Co 



(61) 



The ambiguity of Ci and C2 leads to the gauge modes in 
density and velocity perturbations [28| . 



px a 

-.^ . 

a 



(62) 



where s indicates that the perturbations are confined on 
different Synchronous coordinates. The condition ^ = 
B = does not yield unambiguous results and additional 
definitions are called for. 

Usually, C2 can be obtained by fixing the initial cur- 
vature perturbation |28!], while one gets Ci fixing the 
peculiar velocity of free falling non relativistic species in 
the universe. It is usual to set the peculiar velocity of 
the cold DM to be zero Wc = throughout the expansion 
history, 

Vc = Vc = Ci = {) . 

Hence the condition ip — B — Vc — yields unambiguous 
results. 

In non-interacting case, = is a physical choice be- 
cause it satisfies the Euler equation for cold DM peculiar 
velocity. 

v'^ + nvc = 

However, in the coupled case this point should be care- 
fully investigated. The most generic equation of motion 
for cold DM reads 1231 . 



Pc 



Pc 



(63) 



Compared with non-interacting case, two additional 
terms appear on the RHS of above equation. The first 
term refers to the inertial force density produced by the 
varying rest mass of cold DM in the system and the sec- 
ond term refers to the non-gravitational force density. 
The non-gravitational force consists of two parts. 



SQpc = SQ'p^ \t +Q% 



(64) 
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One is the viscosity Q^Vt due to the momentum transfer 
in DM and DE, where Vt is the average energy transfer 
velocity; another one is external non-gravitational force 
^QpX It due to the collision effect at the early time of the 
universe. If we neglect the non-gravitational force 6Qpc 
in Eq. cold DM particles only suffer the attraction 
of gravity without other external non-gravitational force. 
Eq. still has "free falling" solution Vc = which is 
the same as the non-interacting case. Setting Vc — 0, with 
a completely specified gauge condition, the synchronous 
gauge is as valid as any other gauge. 



The choice of Synchronous gauge has a very close tie 
with the weak equivalence principle. The Synchronous 
coordinate should be chosen to rest upon the local inertial 
frame where the four acceleration A = tp = of observer 
is zero. As discussed above, if there is only gravitational 
forces acting on cold DM bulk, the cold DM particles are 
still "free falling" and the Synchronous gauge is valid for 
the cold DM frame which, in turn, means that cold DM 
frame is a local inertial frame. Since the weak equivalence 
principle is valid in the local inertial frame, it should be 
valid in the cold DM frame. 
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